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FIRST  ORDER  AUTOREGRESSIVE  GAMMA  SEQUENCES 


AND  POINT  PROCESSES 


* 

D . P . Gaver 
and 

* 

P.  A.  W.  Lewis 

Naval  Postgraduate  School 
Monterey,  CA  93940 


1.  INTRODUCTION 

The  Poisson  process  is  a basic  model  for  point  processes 
(series  of  events)  and  can  be  characterized  as  a process  in  which 
tne  intervals  between  events  are  independent  and  exponentially 
distributed.  In  one  of  the  earliest  papers  on  point  processes 
(Wold,  1948)  an  attempt  was  made  to  generalize  the  Poisson  process 
by  obtaining  dependent  but  marginally  exponentially  distributed 
intervals  between  events.  Similarly  Cox  (1955)  attempted  to 
obtain  a sequence  of  random  variables  with  conditionally  exponential 
distributions.  Neither  of  these  attempts  to  generalize  the  Poisson 
process  led  to  analytically  tractable  results.  In  principle, 
of  course,  it  is  simple  to  generate  a sequence  of  marginally 
exponentially  distributed  random  variables  with  Markov  dependence 
if  a bivariate  exponential  random  variable  is  available.  However, 
despite  the  recent  discovery  of  many  bivariate  exponential  random 
5 
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variable  generation  schemes  [Marshall  and  Olkin  (1967),  Hawkes  (1972);  } 

Gaver  (1972);  Downton  (1970)],  none  seem  to  have  simple  enough 
properties,  conditional  and  otherwise,  to  lead  to  analytically 

and  computationally  tractable  models  for  a non-independent  ' 

(Markovian)  and  easily  simulated  sequence  of  marginally  exponentially 
distributed  random  variables.  Interest  in  the  aforementioned 
processes  also  may  have  been  damped  by  the  development  of  physically 
motivated  point  process  models  such  as  the  Poisson  cluster  process 

J 

and  the  self-exciting  process,  for  example  see  Neyman  and  Scott  I 

(1972) . Unfortunately,  many  of  these  point  process  models  are  i 

somewhat  awkward  to  handle  analytically  or  computationally.  j 

In  this  paper  we  show  that  if  one  starts  with  the  usual 
linear,  additive  autoregressive  equation  (first-order  stochastic 

difference  equation)  ] 

J 

Xn  = pXn-l  + en  ' n = 0,  + 1,  + 2,  . . . (1.1) 

J 

where  the  innovation  sequence  (en)  is  one  of  i.i.d.  random  vari-  j 

ables,  then  there  is,  for  0 £ p < 1,  a distribution  for  the  £n’s  j 

■ 

such  that  the  X 's  have,  marginally,  an  exponential  distribution.  J 

n 

The  resulting  exponential  autoregressive  process  (EAR1)  has  several 
attractive  features.  First,  the  sequence  is  obtained  as  an  additive 
random  linear  combination  of  random  variables,  and  is  easy  to 
simulate,  i.e.  to  realize  on  a computer.  In  simulating  a queue, 
for  instance,  one  can  easily  obtain  sequences  of  correlated  service 
times  or  interarrival  times,  with  respectively  an  i.i.d.  exponential 
sequence  or  Poisson  process  as  a special  case  (p  = 0) . These 
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correlated  sequences  are  useful  for  checking  for  the  sensitivity 
of  standard  queueing  results  to  departures  from  independence,  and 
for  other  robustness  studies. 

A second  feature  is  that  the  EARl  process  (1.1)  is  analyti- 
cally tractable;  one  can,  for  instance,  obtain  the  Laplace  transform 
of  the  distribution  of  sums  of  the  random  variables.  This  is 
essential  if  the  process  is  to  be  used  as  a model  for  point  pro- 
cesses and  one  is  to  obtain  the  point  spectrum. 

A third  feature  of  the  EARl  process  is  that  its  structure 
leads  to  an  extended  model  for  a sequence  of  marginally  exponentially 
distributed  random  variables  with  the  correlation  structure  of  a 
mixed  autoregressive  moving  average  process  (ARMA  p,q) , the  orders 
of  the  autoregression  and  the  moving  average  being  p and  q 
respectively.  Various  properties  of  this  exponential  process, 
called  the  EARMA(p,q)  process,  are  detailed  in  Lawrance  and  Lewis 
(1977) , Jacobs  and  Lewis  (1977)  and  Lawrance  and  Lewis  (1978) . 

In  Section  2 of  this  paper  we  discuss  the  exponential 
solution  of  equation  (1.1),  as  well  as  questions  of  stationarity 
and  mixing.  The  distribution  of  sums  of  random  variables  from 
this  process  is  discussed  next;  this  relates  in  particular  to  the 
use  of  the  process  to  model  intervals  between  events  in  point 
processes.  Joint  distributions  and  conditional  correlations  of 
two  and  three  random  variables  are  derived  in  Section  4 . The 
estimation  of  the  exponential  parameter  1,  the  reciprocal  of 
E(Xn),  and  the  correlation  parameter  p are  briefly  considered 
in  Section  5 for  a fixed  or  random  number  of  observed  random 
variables.  It  is  shown  that,  because  of  a certain  degeneracy 
in  th€'  process,  p can  be  estimated  exactly  in  a long  enough  sequence. 
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Other  solutions  to  equation  (1.1)  are  considered  in  Section  6. 


In  fact  random  variables  for  which  solutions  to  (1.1)  (in  trans- 
formed version)  exist  are  known  as  self -decomposable  random 
variables  or  random  variables  of  type  L;  see  Feller  (1971,  p. 
588-90) . Two  cases  of  interest  in  modelling  point  processes 
are  discussed.  These  require  consideration  of  Gamma-distributed 
intervals,  and  intervals  with  a mixed  exponential  distribution; 
questions  connected  with  the  latter  remain  unresolved. 

Finally  in  Section  7 we  consider  the  question  of  obtain- 
ing autoregressive  and  Markovian  exponential  sequences  with 
negative  correlation.  It  may  be  seen  that  (1.1)  has  no  solution 
if  P is  negative.  However,  as  soon  as  multivariate  exponential 
sequences  are  considered — in  particular,  multivariate 
sequences  which  are  antithetic  realizations  of  each  other 
— we  discover  a method  for  obtaining  processes  with  exponential 
or  other  specified  marginals  and  that  exhibit  negative 
correlation. 
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2.  EXPONENTIAL  SEQUENCES 

In  what  follows  {E  } will  always  be  a seouence  of  i.i.d. 

n 

exponentially  distributed  random  variables  with  parameter  X? 


P { E < x)  = F„ 
n — 


(x)  = 1 - e 


n 


-Xx 


(x  > 0,  X>  0) 


(2.1) 


= 0 


(x  < 0) 


Also  an  exponential  (A)  random  variable  will  mean  a random 
variable  with  distribution  (2.1). 


2.1.  The  exponential  first-order  autoregressive  sequence  (EAR1 ) . 

The  starting  point  of  the  work  in  this  paper  is  the  question 
as  to  whether  the  autoregressive  equation  (1.7) 


X = 
n 


pX  , + e 
K n-1  n 


= I 

j=0 


p3e 


n-3 


Ipl  < 1 


(2.2) 


has  a solution  for  a given  distribution  of  X^  (note  that  the 

expansion  in  (2.2)  as  an  infinite  moving  average  is  valid  because 

|p|  < 1).  Now  xn_i  is  a function  only  of  en_i'  £n-2'  **'  an£* 

is  therefore  independent  of  en>  Therefore  the  Laplace-Stielt jes 

transform  4>v  (s)  of  the  distribution  of  X is,  at  least  for 
n n 

s > 0 , 


4>x  (s)  = E[exp(-sXn)]  = E [exp(-spXn_1  + en>  ] 


= 4>v  (ps)  4)  (s)  . 

n-l  n 
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Thus 


<J>  (s)  = 


(S) 

n 


n 


x 


(ps) 


(2.3) 


Assuming  that  the  Xn  sequence  is  marginally  stationary,  we  get 


the  basic  equation 


0E(S) 


X 


(s) 


<f>x  ( ps) 


Now  it  is  clear  that  if  we  require  the  Xn  to  be  positive  random 
variables,  then  if  p is  negative  so  is  pXn  and  we  need  the 
error  term  which  is  independent  of  xn_]/  to  maXe  xn  positive. 

Thus  for  positive  random  variables  there  will  clearly  be  no  solution 
to  (2.3)  for  p < 0.  However  if  we  let  0 <_  p < 1 and,  side- 
stepping the  general  question  of  existence  and  uniqueness  of 
solutions  to  (2.3),  require  the  xn's  to  be  exponential  with 


Vs)  - A-hr  ' 


then  equation  (2.3)  forces 

A \ X X + ps 

Vs)  = m • — — 


X + ps 
X + s 


(0  < p < 1,  X > 0) 


= p + (1-p) 


X + s 


(2.4) 


The  latter  expression  is  the  Laplace-Stielt jes  transform  of  a 

random  variable,  in  fact  a non-negative  random  variable  which 

has  an  atom  of  mass  p at  zero  and  which  is  exponential ( X)  if 

positive.  Thus  we  can  write  the  difference  equation  generating 

the  series  {X„}  as 
n 


4 


; 
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X = pX  , + e 
n n-1  n 


| Pxn_i  W*P*  P 

l Pxn_!  + En  w.p.  (1_P) 


0 < p < 1 


(2.5) 


= pXn-l  + InEn' 


(2.6) 


where  {1^}  is  an  i.i.d.  sequence  in  which  In  = 1 with  prob- 
ability ( 1— p ) , I = 0 with  probability  p and  {En}  is  an 
i.i.d.  sequence  of  (X)  exponentials.  It  is  easily  verified, 
directly  from  (2.2)  or  from  the  definition  of  e , that 


EUnl  = (1-P)  x 


var(e  ) = \ (1-p2) 
n > ^ 


(2.7) 


There  are  several  points  to  be  made  about  the  sequence  Xn: 

(i)  When  p = 0,  then  {Xn>  is  an  i«i»d.  sequence  of 
exponential (X)  r.v.s. 

(ii)  The  representations  (2.5)  and  (2.6)  of  the  process  as 

an  (additive)  random  linear  combination  of  X , and  E , 

n-1  n 

which  has  a distribution  independent  of  p,  are  sometimes 
more  convenient  than  is  (2.2). 

(iii)  Although  Xn  is  strictly  a linear  process  (AR1 ) , it 

differs  from  the  normal  or  Gaussian  process  ((2.2)  with 
en  normally  distributed  of  Gaussian)  in  that  the  mean , 
variance  and  higher  moments  of  the  en ’ s are  functions 


7 


(iv) 


of  the  parameter  p.  Thus  one  cannot  apply  general  theorems  about 

linear  processes  to  the  exponential  sequence  because  such 

theorems  usually  assume  that  the  e ’ s have  a distribution 

n 

which  is  free  of  the  parameter  p.  As  an  example,  Xn  is 

exponentially  distributed  even  when  p is  close  to  one, 

despite  theorems  for  linear  processes  which  assert  that 

in  this  case  the  are  approximately  normally  distributed. 

For  instance,  if  the  en's  are  independent  of  p then 

the  skewness  of  the  stationary  distribution  of  X^  can 

be  shown  to  approach  zero  as  p + 1;  in  the  special  case 

of  (2.4),  this  quantity  is  independent  of  p and  equal 

to  2--the  value  associated  with  the  exponential  d.f. 

There  are  many  other  sequences  of  random  variables  {X^} 

with  exponentially  distributed  marginals  and  the  first-order 

Markov  property.  In  fact,  given  any  bivariate  exponential 

distribution  F„  _ (x,,x„),  assumed  for  convenience  to 
E1,E2  1 2 

be  absolutely  continuous,  with  conditional  density 

f£  |E  (x2;x^),  then  we  construct  a Markovian  sequence  with 


joint  density  for,  say,  Xn,  Xn_^, 


, XQ  as 


f\ V*0 Xq) 

fE  |E  (xn;  Xn-1)  f E |E  (xn-l;  Xn-2} 

n n-1  n-1  n-2 


fE1lE0(Xl,3C°)  fE0  (Xq)  ’ 
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An  interesting  and  useful  feature  of  the  EARl  sequence 
given  by  the  solution  (2.6)  of  the  autoregressive  equation  is 
that  it  is  a (random)  linear  combination  of  i.i.d.  exponential 
sequences,  and  is  thus  easy  to  simulate  on  a computer,  gives 
reasonably  tractable  analytic  results  and  in  turn  suggests  the 
definition  of  processes  with  more  complicated  correlation 
structure  and  exponential  marginals  (Lawrance  and  Lewis,  1977, 
Jacobs  and  Lewis,  1977,  Lawrance  and  Lewis,  1978) . The  correlation 
structures  are  essentially  the  same  as  those  for  linear  processes 
(Lawrance  and  Lewis,  1978). 


2.2.  Serial  correlations,  stationarity  and  mixing 

Simple  computations  show  that,  as  is  true  of  any  regular 
Markov  process,  the  serial  correlations  for  the  EARl  process  are 


c.  = corr(XnXn+j)  - pj  , 


0 < p < 1. 


Note  that  the  correlations  are  always  positive.  The  spectrum  of 
the  sequence  is 


f+(u> 


1 

)=  — {1  + 2 l c.  COS  ( jco  ) } ( 0 < U)  < TT  ; 0<p< 

j=l  3 - - - 


1) 


1 

IT 


1 - p- 


1 + p - 2p  cos 


U) 


For  p = 0 we  have  f+(w)  = 1/tt  since  the  sequence  is  i.i.d. 


Moreover  it  is  not  difficult  to  show  that  the  process  is 


(2.8) 

(2.9) 

with  e defined  by  (2.4).  This  is  no  surprise,  for  the 
n 

process  is  Markovian  and  is  constructed  to  have  an  exponential 
distribution  in  the  stationary  case. 

With  Eq  = Xg  it  is  shown  in  Jacobs  and  Lewis  (1977)  that 
the  sequence  {Xn>  is  strong  mixing  in  the  sense  of  Rosenblatt 
(1971)  . 


stationary  if  one  takes 


X0  " E0 


X = pX  . + e = pX  . + I E 
n H n-1  n H n-1  n n 


3.  SUMS  OF  STATIONARY  INTERVALS 

An  important  aspect  of  the  stationary  sequence  {X^}, 

especially  when  it  is  used  to  model  the  intervals  in  a point 

process,  is  the  moments  and  distributions  of  sums 

T = X +X  ,+•••+  X ..  The  Laplace-Stielt jes  transform 
r n n+1  n+r-1 

of  T^  is  found  directly  by  noting  that,  in  terms  of  Xn, 
(2.2)  may  be  modified  to  give  this  representation: 
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It  follows  that 


(3.2) 


If  the  process  starts  with  the  stationary  distribution,  i.e., 
if  is  exponential ( A) , then 


Clearly  the  distribution  of  can  be  found  explicitly  by 

inverting  (3.3),  which  can  be  accomplished  by  expanding  in 

partial  fractions.  So  also  can  an  expression  for  the  interval 

spectrum.  The  result  is  analytically  awkward,  and  will  not 

be  quoted.  In  order  to  obtain  (i)  the  intensity  function,  m^(t), 

of  a point  process  with  EARl  intervals,  (ii)  the  point  spectrum 

g+(u),  and  (iii)  the  distribution  of  the  counting  process  N (t) 

(for  definitions  see  Cox  and  Lewis,  1966,  Sec.  4.5)  one  must  be 

able  to  sum  $ (s)  over  r;  there  appears  to  be  no  neat 

r 

explicit  formula  for  any  of  these  functions.  The  variance  of 
Tr  may  be  calculated  directly  from  (3.1)  since  the  innovations 
{ e n } are  independent: 


Hence  for  p < 1 the  index  of  dispersion  (Cox  and  Lewis  (1966) , 
p.  71)  is 


4.  JOINT  DISTRIBUTIONS  FOR  THE  STATIONARY  SEQUENCE 

Any  pair  Xn,  Xn+r  in  the  stationary  EARl  sequence  has 
a bivariate  exponential  distribution.  Consider  and  xn+^. 

Directly  from  the  definition  (2.5)  and  (2.7)  it  can  be  shown 
that  the  conditional  random  variable  xn+^»  given  Xn  = x,  has 
an  atom  of  mass  p at  px;  otherwise  it  is,  with  probability  (1-p) , 
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an  exponential ( X)  random  variable  shifted  px  from  the  origin. 
The  regression  of  xn+i  on  Xn  = x 


E(Xn+l'Xn  - x) 


x + (1-p)  J , 


(4.1) 


2 2 

which  is  linear.  Moreover  var(Xn+1IXn  = x)  = (1-p  )/X  , a 

constant  independent  of  x,  and  xn+1  is  never  smaller  than 

pXfi.  The  Laplace-Stielt jes  transform  of  the  bivariate  exponential 

distribution  of  X Ll  and  X is 

n+l  n 


;X 


n+1 


S2> 


- E{exP(~SlXn  - s2Xn+1)> 

= ECexpf-s^^  + p s 2 ) XR ] } <Pe  (s2) 

X ( . . X | 

+ s,  + ps„  )p  ' p X + s,,  ( 


X + s^  + p&2 


p X + s1  + ps2  + ( X + s^  + ps2)(x  + s2  ) 


(4.2) 


The  bivariate  distribution  of  {Xn+-^,  X^}  has  a singular  component 

along  the  line  XR+^  = Xn»  with  probability  p (the  first  term 

in  (4.2),  and  a continuous  component  in  the  space  X^+^  > pX^ 

(the  second  term  in  (4.2))  which  is  the  joint  distribution  of  Xfi 

and  px„  + E . Note  that  the  distribution  is  not  symmetric  in 
n n 

X . and  X , since  the  transform  <J>V  v (s.  ,s~)  is  not 
n+l  n xn+l ' xn  1 * 

symmetric  in  s^  and  s2«  This  simply  means  that  the  process 
is  not  time-reversible,  as  is  the  normal  or  Gaussian  ARl  process. 

This  asymmetry  appears  in  the  conditional  mean  of  Xn, 
given  Xn+1  * x,  which  is  obtained  from  (4.2),  as 

This  asymmetry  also  becomes  evident  in  the  non-linear  form 
of  the  conditional  mean  of  Xfl,  given  xn+l  - x'  which  is  obtained 
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from  (4.2),  as 


- m^pT  {1  - exP[-  <?  - 11  X*1  > * E>VlIXn-x)  • 

The  conditional  variance  of  X , given  X , . = x,  can  also  be 

n n+1 

obtained;  unlike  var(Xn+1IXn  = x)  it  is  not  a constant. 

Higher  order  directional  moments  of  xn+^»  xn  are  given 
in  Jacobs  and  Lewis  (1977) , and  again  evidence  the  directionality 
of  the  process: 


2,1(1) 

- E(x^ 

W 

- E'Xn' 

E(Xn+l> 

1,2^) 

= E (X 

n 

- E(Xn> 

E(xn+1> 

. 2p  (1+p) 
— 

Since  the  process  is,  from  its  definition,  Markovian,  the 

conditional  correlation  of  X and  X„  , , given  X = x,  is 

n+i  n-i  n 

zero.  We  do  not  discuss  distributions  of  triples  any  further. 

In  principle  it  is  simple  to  write  down  transforms  of  the  joint 
distribution  of  any  set  of  k xn's*  one  obvious  use  of  this 
result  would  be  in  obtaining  the  distribution  of  the  sum  of 
Xn,  ...  , xn+jc_2 • However,  this  was  obtained  by  a direct  argument 
in  the  previous  section. 


5.  ESTIMATION  AND  DEGENERACY 

Although  it  is  not  possible  to  write  down  the  likelihood 
equation  for  the  EAR1  process  in  a tractable  form  where,  say,  we 
observe  X^,  ...  , X^j,  it  seems  at  first  sight  that  good  estimates 
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T 


of  1/ A and  p can  be  obtained  from  the  sample  mean,  X = T^/n, 
whose  variance  is  given  at  (3.4),  and  a first-order  serial 
correlation  coefficient  estimate  (3^  with  modifications  to  take 
account  of  the  fact  that  the  variance  equals  the  mean  squared, 
i.e.  Var(X)  = l/A^. 

However,  there  is  a degeneracy  in  the  process  which  makes 
it  possible  to  find  p exactly  in  a long  enough  run  of  xn's* 
and  to  estimate  A more  precisely  than  is  possible  using 
straightforward  moment  methods. 

To  identify  the  degeneracy,  which  we  call  the  zero-defect, 
we  note  that  in  the  process  there  are  runs  of  xn's  which  are 
equal  to  the  previous  value,  xn_^ f times  p.  Moreover,  given 
that  there  is  a run  of  length  R of  this  type,  R has  a "geometric 
plus  one"  distribution  with  parameter  p; 


] 

i 


P(R  = i>  = (1  - p)p1_1  , i = 1,2,3, 


(5.1) 


E{R)  - ? 


var (R) 


This  behavior,  or  a tendency  towards  it  in  data  would  be  very 

evident  in  plots  of  a sample  sequence  of  Xr  values. 

Thus  to  estimate  p in  an  observed  series,  let  Z = X ,/X  , 

n n+1'  n ' 

for  n = 2,3,  ...  . Then 


= P 


if  Xn+1  = pXn  W‘P*  p 

if  Xn+1  f pXn  w.p.  (1-p)  . 
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The  probability  that  ZR  = p,  or  Zn  / p,  is  independent  of  previous 

Zn' s and  P(Zn  = p)  = p.  Hence  if  we  observe  the  minimum  value  of 

Zn,  n = 2,3,...  until  we  get  a tie,  that  tie  value  is  p.  The 

time  to  obtain  this  tie  is  the  sum  of  two  "geometric  plus  one" 

random  variables,  with  total  mean  2/p  and  total  variance 
2 

2(l-p)/p  . The  implication  is  that  we  can  find  p exactly  after 

a random  number  of  observations  for  the  present  model.  Moreover 

since  one  then  knows  those  X ’ s which  are  made  up  of  px  . 

n n-l 

plus  an  En,  it  is  possible  to  compute  the  En's  exactly  (except 
for  two)  and  estimate  E(X)  = 1/X  directly  from  the  observed  En's* 
For  a fixed  length  observation  of  the  process,  say 
X^,  ...  , X^ , the  probability  that  p can  be  obtained  exactly 
is  the  probability  that  the  sum  of  two  "geometric  plus  one" 
random  variables  is  less  than  or  equal  to  N.  If  two  minimum 
values  are  not  observed  by  time  N,  a good  estimate  of  p will 
be  the  minimum  of  the  Z^'s.  This  will  either  be  p,  or  p plus 
a small  bias  term.  The  mean  can  be  estimated  as  the  sample 
average  in  the  usual  way. 

This  degeneracy  in  sequences  of  positive  random  variables 
generated  from  the  stochastic  difference  equation  seems  to  be 
inherent  in  the  present  procedure.  Slightly  more  realistic  but 
complicated  first-order  autoregressive  exponential  processes 
without  this  zero  defect  will  be  discussed  elsewhere. 
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6.  OTHER  SOLUTIONS  TO  THE  AUTOREGRESSIVE  EQUATION 

Random  variables  X for  which  the  transformed  first-order 
autoregressive  equation 


<t>x(s) 

~r ? — — r = <f>  (s) 

<t>x(ps) 


(6.1) 


has  a solution  <)>e  (s)  for  each  0 < p < 1 which  is  the  transform 
of  a distribution  function  are  called  random  variables  of  class  L 
(Feller,  1971,  p.  588),  or  self -decomposable  random  variables. 
Although  there  has  been  much  recent  interest  in  finding  class  L 
random  variables,  the  connection  with  the  first-order  autoregressive 
process  does  not  seem  to  have  been  made.  On  the  other  hand  there 
have  been  some  attempts  to  find  non-normal  solutions  to  the  auto- 
regressive equation  without  connecting  them  to  class  L theory 
(see,  e.g.  Bernier,  1970) . 

The  limitation  of  the  theory  of  class  L random  variables  as  it 
relates  to  the  present  work  is  that  it  requires  a solution  for  each 
0 < p < 1,  which  is  the  case  for  exponential ( A ) random  variables  X. 
This  full  range  of  p is  desirable,  but  may  not  occur. 

We  do  not  explore  the  full  connection  here,  but  consider  only 
two  types  of  random  variables.  Gamma  distributed  random 
variables  and  mixed  exponential  random  variables.  These  are 
frequently  used  as  alternatives  to  the  exponential  in  modelling 
stochastic  phenomena  such  as  response  times  in  queues. 
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6.1.  The  Gamma  Autoregressive  Process  GARl . 


A necessary  condition  for  X to  be  in  the  class  L is 
that  it  be  infinitely  divisible,  and  therefore  we  consider 
Gamma  (A,k)  variables  with  density 

,k  k-1o- 

f Y ( x ) = - (A  > 0;  k > 0;  x > 0)  (6.2) 

x kkr(k) 

and  Laplace-Stielt jes  transform  (of  the  stationary  distribution) 

‘6-3> 

Then  the  solution  to  (6.1)  is 

- ( *-4ir)k  * (»  + rrr)k  (6-4) 

It  follows  from  a result  reported  by  Feller  ((1971),  p.  452)  or 
directly  from  Theorem  1,  p.  450)  that  [(A  + ps)/(A  + s) ] is 
the  transform  of  an  infinitely  divisible  distribution,  and  hence 
(6.3)  is  actually  the  transform  of  (an  infinitely  divisible)  distri- 
bution function  for  every  real  k > 0.  Thus  we  can  in  principle 
generate  an  autoregressive  process  with  gamma  marginals  (6.2) 
by  utilizing  the  {en)  process  characterized  by  (6.3).  Here 
are  a few  simple  special  cases. 

k = 2:  ♦E(s)  = p2  + 2p(i-p)  j-i-x + a~p>2  (rrrf1  (6-5) 

k = 3=  *EU>  - P2  + 2 p2  ( 1“P)  j-fj  + 2 p ( 1-p)  3 ( 5-t-t:) 

♦ «-»> 3 ( rfr)3  • (6'6) 
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Thus  for  k = 2 the  random  variable  e may  be  thought  of  and 
realized  as  a convex  mixture  of  a degenerate  random  variable 
with  mass  at  zero  and  an  exponential ( X)  and  a Gamma  (X,2) 
distribution.  Alternatively  if  and  E ^ are  independent 

exponential (X)  random  variables,  the  sum  is  chosen  with  probability 

2 2 
(1-p)  , E^  is  chosen  with  probability  2p(l-p)  , and  neither  of 

2 

E.^  or  E2 , but  rather  zero,  is  chosen  with  probability  p . 

Unfortunately,  no  really  simple  way  of  generating  random  vari- 
ables with  Laplace-Stielt jes  transform  (6.4)  for  noninteger  k is 


known,  though  representation  as  an  infinite  sum  of  weighted 
Gamma  variables  is  possible  (Bernier  (1970)).  Moreover  for  a 
case  of  much  interest  k < 1,  so  that  X is  overdispersed 
relative  to  an  exponential  random  variable,  and  the  degeneracy 
or  zero-defect  of  weight  p becomes  very  prominent. 

As  was  true  for  the  exponential  case,  we  have  for  the 
correlations,  which  are  all  positive, 

Cj  = p^  0 <_  P < 1; 

also 


E (e)  = X(l-p) 

Var(e)  = X2(l-p2)/k  . 


As  k -*•  00 , the  degeneracy  disappears  and  X tends  to  become 
normally  distributed. 
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6.2.  The  Mixed  Exponential  Process  MEARl 

Consider  now  X to  be  a mixture  of  two  exponential  densities 

-Ax  -A2x 

f (x)  = TT^A-^e  + Tr2^2e  ' 

x > 0,  = 1 — tt2  0,  A^  < A 2 . (6.7) 


For  = 0 , ir^  = 1,  and/or  A^  = A2  this  of  course  reduces 

to  an  exponential  density.  There  are  two  main  cases  to  consider, 
besides  these  special  cases. 

If  0 <_  tt ^ <_  1,  then  we  always  have  a convex  mixture  of 
exponential  densities:  fx(x)  is  a ProPer  density  function  and 
the  coefficient  of  variation  of  X lies  between  1 and  infinity 
(Cox,  1964) , so  that  X is  overdispersed  relative  to  an 
exponential  random  variable;  X may  be  generated  as  a mixture 
of  two  exponential  random  variables  with  parameters  A^  and 
A2,  and  X is  infinitely  divisible  (Feller  (1971),  p.  452). 

If  tt^  is  greater  than  1,  so  that  tt2  is  negative,  a 
necessary  and  sufficient  condition  for  f^(x)  to  be  a p.d.f. 
is  that  + ^2^2  — 0 or  ec3uivalentlY  that  tt1  <_  [1  - (Aj/A2)]  1 . 

Here  is  alternative  way  of  writing  the  transform  of  X. 
First,  and  directly. 


♦XU)  = 


X1 + 


+ IT 


2 A-  + s 


(6.8) 


Now  since  factor  to  obtain 
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(s)  = 


^1  / X2  + 5 \j 

"1  X2  \ Xx  + s 

j1  - "1  l1  - Tj)  + ”l( 


1 ~ T~ 


V X1  + S 


(6.9) 


Now  even  if  tt  > 1 (tt  < 0)  but 
1 2 


t ■ "i  l1  ■ ^) 


< i , 


(6.10) 


or  equivalently 


7T2X2  + 7TiXl  X 0 


(6.11) 


then  ( s ) is  the  transform  of  the  sum  of  an  exponential  (X 

Jrk. 

random  variable  and  a random  variable  having  the  distribution 
of  the  e-innovation  of  an  EARl  process,  see  (2.4).  It  follows 
that  X is  infinitely  divisible. 
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. It  is  thus  worth  investigating  whether  one  can  generate 
an  autoregressive  process  (Xn)  with  mixed  exponential 
marginals  for  some  or  all  p,  the  latter  being  the  question 
as  to  whether  the  mixed  exponential  random  variable  is  in  the 
class  L.  We  have,  directly 


<Pe  (s) 


V5) 

4>x  (ps) 


*^1^1  ^2^2 
s +A^  + s +A  ^ 

tt  A i ^2^2 
ps  + A^  + ps  + A2 


( A ^ + ps)  ( A 2 + ps)  (A1A2  + tt-^A^s  + 7T2'^2S^ 
( A ^ + s)(A2  + s ) ( A 1 A 2 + 7T1A1ps  + tt  2 A 2 p s ) 


After  considerable  simplification  this  Laplace-Stielt jes 
transform  can  be  written  as  a mixture  of  a degenerate  random 
variable  with  mass  p at  zero,  and  with  probability  (1-p) 
(possibly)  a random  variable  Y which  has  (possibly)  density 
function 


-A,x  -A„x  . 

fy(x)  = Vle  + Y2X2e  + Y3  ^ e xc  p 


(6.13) 


where  y^  + y2  + Y3  = 1,  A^  < > 0,  1 and 


C = 


*1*2 

^l^l  + <1_1ri)A2 
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U2  - pxx)  (C  - Xx) 


(6.14) 


Y1  = 


(X2  - X1) (C-  pXx) 


= PX  ; 


Y2  " 


(X2  - C)  (pX2  - Xx) 


(X2  - Xx)  (pX2  - C) 


= P2 


(X1  -pX2) 
(C  - pX,) 


•;  (C  ^ pX2) 


= 


(C  ' V (X2  ' C)  1 1 

fc  -”xjl  X (C  - pX2r  - P3  X (C  - pX2)  ; 


(6.15) 


(6.16) 


Now  if  is  restricted  to  be  0 < < 1,  then  X^  < C < X2 


and  p^  > 

0,  p 

2 > 

0,  p3  > 0.  Then  it  can  be  shown  that  the 

signs  of 

V 

CM 

Y3  depend  on  whether  0 < p < X^/X^, 

x.j/x2  < 

P < C/X2 

or  C/X2  < p < 1: 

(i) 

Yi 

is 

always  positive; 

(ii) 

Y2 

and  are  positive  if  0 < p <_  X^/X2; 

(iii) 

Y2 

is 

negative,  y3  is  positive  if  ^ l/^ 2 < P < C/X^; 

(iv) 

^2 

is 

positive,  y3  is  negative  if  C/X2  < p < 1. 

While  this  result  is  relatively  simple,  the  only  case  in  which  it 
has  been  possible  to  establish  that  fY<x)  is  a p.d.f.  is  that 
in  (ii)  , i.e.  0 < p X^/Xj.  The  ratio  X^/Xj  is  rou9hlY 
related  to  the  overdispersion  of  X relative  to  the  exponential 
distribution}  the  smaller  the  value  Xj,/X2,  the  greater  the 
dispersion  and  the  smaller  the  admissible  range  of  p.  It  seems 
probable  that  fY(x)  is  not  a density  function  for  all  P;  at 
all  events  efforts  to  prove  that  it  is  in  the  class  L using 
(6.10)  and  characterization  theorems  (Feller,  1971)  have  failed 
at  the  time  of  this  writing.  A complete  understanding  is  not  now 
available. 
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6.3.  Discussion 

Note  the  differences  between  the  mixed  exponential  auto- 
regressive sequence  and  the  Gamma  (k  < 1)  autoregressive  sequence; 
the  former  appears  to  require  a restricted  range  for  p,  e can  be 
simulated,  and  the  zero  defect  has  probability  P;  the  latter 
sequence  is  valid  for  0 £ p < 1,  e cannot  at  present  be 
simulated,  and  the  zero-defect  has  probability  p > p . 

The  magnitude  of  this  zero-defect  alone  probably  makes  the 
Gamma  process  less  useful  than  the  mixed  exponential  process. 

The  tail  behavior  of  the  marginal  distributions  is  also 
different. 

Laplace- Stie It jes  transforms  of  distribution  of  sums, 

T of  X ' s are  obtained  from  (3.2)  for  both  the  Gamma  and 
rl  n 

mixed  exponential  autoregressive  sequences.  These  cases  are 
still  under  study. 
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7. 


MULTIVARIATE  SEQUENCES  AND  NEGATIVE  CORRELATION 


T 


Even  though  the  autoregressive  equation  (1.1)  has  no  solution 

representing  a positive  random  sequence  when  p < 0 it  is  possible 

to  modify  it  so  as  to  produce  negatively  correlated  exponential 

sequences.  These  will  not  be  individually  Markovian,  but  will  have 

Markovian  properties  in  a bivariate  sense.  A suggestion  as  to  how 

to  proceed  comes  from  writing  Xn  in  terms  of  the  innovations 

(e„ } as  in  (2.2): 
n 


Xn  = 


n 

3I0 


p3e 


n-j 


p£n-l  + 


2 

P e 


n-2 


+ • 


pn£0 


(7.1) 


Now  consider  generating  two  (or  more)  series  in  parallel,  using 
independent  pairs  of  innovations,  (£n»  e^) # for  which  if  desired 
Cor(en,  e^)  < 0;  e.g.  for  0 £ p < 1 and  starting  with  the  first 
innovation. 


Model  1:  X = pX' 

n n-1 

+ e = 
n 

e + 
n 

peA-l 

+ 

p2£n-2 

+ 

p3e ' , 
p n-3 

+ • • • 

XA  - pXn-l 

+ eA  = 

e 1 + 
n 

pen-l 

+ 

p2pA-2 

+ 

P\-3 

+ • • • 

Model  2:  X 

n 


= PX, 


n-1 

+ 

en  = 

en 

+ pen-l 

+ 

2 

p en-2 

+ 

3 

P e ■, 

n- J 

+ • * • 

’ 1 

n-1 

+ 

e ' = 
n 

eA 

+ p£i-l 

+ 

p^E  ' 

p en-2 

+ 

°3eA-3 

+ • * • 

(7.30 
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In  order  for  {Xfi}  and  {X^}  to  exhibit  asymptotically 

exponential  marginal  distribution,  it  is  clear  that  the  marginal 

innovation  sequences  (e  } and  {e'}  must  themselves  be  iid 

n n 

with  an  atom  p at  zero,  and  otherwise  an  exponential  distribution. 
This  can  be  accomplished  in  several  ways;  here  are  two: 


where  {Un>  is  a sequence  of  uniform  (0,1)  "random  numbers"; 

(E  } and  {E'}  are  then  called  antithetic,  see  Hammersley 
n n 1 

and  Handscomb  (1964)/  and  are  maximally  negatively  correlated, 

2 

having  correlation  1 - (it  /6)  = - 0.6449;  see  Moran  (1967). 
Each  series  in  Models  1,  2 receive  independent  innovations  of 
exactly  the  type  described  in  Section  2 and  that  were  there 
shown  to  lead  to  exponential  marginals  as  n °°.  In  this 


E[e,e']  = (1-p)  -i-  / £n  u £,n(l-u)du  = (1-p)  -i-  [* 2 - 

ao  v L 


2 

Cov ( e , e ' ) = (1-p)  “T  [2  - V]  “ ~ 


Also  Cov(e  ,c’)  >0  if  p > (n  /6  - 1) . 
n n 


(7.5) 

(7.6) 
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then  in  true 


(ii)  If  the  distribution  function  of  e is  F , 

n en 

antithetic  fashion  determine 


e = F-1(U  ) ; 
n e n 


In  particular,  if 


£'  = F 1 ( 1— U ) 
n e n 


(7.7) 


(° 

if 

X 

< 

0 

n 

I 0 

if 

X 

= 

0 

n I 

- Ax 

k P + (1-P) (l-e  AX) 

if 

X 

> 

0 

Then  innovation  leading  to  marginal  exponentials (A)  are 
obtained;  see  Section  2.  In  this  case 


£ =0 
n 

if 

Un  < P*  eA  " 0 

if 

un  > 1-P 

ii 

i 

Hm 

3 

1-°n) 

1-P  / 

if 

Un  > P* 

x £” 

ft)  “ 

un  < 1-P 

(7.9) 

Consequently 

E[£,£;] 

1 

" A2 

i-p 
/ in 
0 

(*)  {&) 

du  , 

0 < p < | 

(7.10) 

= 0 , i-  < p £ 1 


and  hence 
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1 


2 


0 < P i I 

(7.11) 


Scheme  (ii)  is  capable  of  generating  negative  correlations 

of  greater  magnitude  than  is  Scheme  (i) , as  is  clear  from 

examining  the  situation  1/2  < p. 

The  expansions  (7.2)  and  (7.3)  allow  the  computation  of 

the  lagged  cross  covariances  between  X and  X' . For 
” n n 

n -*■  oo  we  find 


Model  1:  Cov (X  ,.,X’) 
n+j  n 

1-p 

Cov(e,e‘)  , 

j = 0,1,2,... 

(7.12) 

-j 

Cov ( e , e ' ) 

for  j = 0,2,4,... 

Model  2:  Cov (X^ ^ . , X* ) 
n+ j n 

~ J 

X 

(7.13) 

1-P 

Var  (e) 

for  j = 1,3,5,... 

Of  interest  is  the  fact  that  the  series  of  Model  2 may  in 
fact  exhibit  negative  correlations;  for,  directly  from  (7.2), 


and  as  n •+• 

- 

Cov(xA+i'xA) 

~ P 

Cov(e, e' ) 

, 2 
1-p 

Cov(Xn+2'V  - 

Cov<XA+2'XA) 

\T  3 y*  ( c~  \ 

1-P2 

var v t ; 

Cov(Vj-V  - Cov(xA+j'xA)  ~^-5T 


Cov ( e , e ’ ) , j 
Var(e),  j 


1.3.5.. .. 

2.4.6.. .. 
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Consequently  if  scheme  (i)  or  (ii)  above,  or  a multitude  of  other 
possibilities,  generate  the  innovation  pairs,  then  the  odd-numbered 
covariances  will  in  fact  be  negative. 

Clearly  the  above  generation  scheme  can  be  generalized, 
e.g.  by  including  more  equations,  {Xn},  {X^},  ^n^  ' etc*  ' anc*  by 
allowing  the  innovations  to  have  different  marginal  properties. 
Further  investigations  are  planned. 
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8.  SUMMARY  AND  CONCLUSIONS 

We  have  presented  a simple,  autoregressive  Markovian  sequence 
{Xn>  of  exponential  variates  EARl  which  is  an  additive  random 
linear  combination  of  the  previous  value,  X and  an  independent 
exponential  random  variable.  The  simplicity  of  this  structure 
allows  one  to  model  in  an  intuitive  way  dependencies  in  stochastic 
systems.  Jacobs  (1978)  has  considered  cyclic  queues  with  EARl 
service  times  and  found  that  the  correlation  may  produce  a 
significant  effect;  more  general  queueing  schemes  which  generate 
multivariate  exponential  sequences  are  given  by  Lewis  and  Shedler 
(1978)  . 

Maxima  of  the  Xn  in  the  EARl  process  have  been  studied 
by  Chernick  (1977) . 

Two  other  marginal  distributions  for  the  X 's  have  been 

n 

studied,  the  Gamma  distribution  and  the  mixed  exponential.  The 
former  is  known  to  be  a type  L,  i.e.  satisfies  ( 6.1)  for  all 
p,  while  the  mixed  exponential  appears  to  satisfy  (6.1)  only  for 
a limited  range  of  p.  Other  type  L distributions  will  be 
investigated  in  the  context  of  the  modelling  of  independent 
stochastic  sequences  elsewhere. 

Extensions  of . the  first-order  autoregressive  structure 
for  exponential  marginals  to  higher-order  autoregressions, 
moving  averages  and  mixed  autoregressive-moving  average  structures 
has  been  given  by  Lawrance  and  Lewis  (1977) , Jacobs  and  Lewis 
(1977)  and  Lawrance  and  Lewis  (1978) . The  possibility  of  extending 
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the  EARl  structure  depends  on  the  fact  that  the  {en>  error 
sequence  is  a mixture  of  a random  variable  with  mass  at  zero 
and  an  exponential  distribution. 

Three  other  possibilities  will  be  detailed  elsewhere — 
a slight  complication  which  brings  in  another  parameter  but  gets 
rid  of  the  "zero-deficiency"  in  the  EARl  process;  extension  of 
mixed  correlation  structures  to  give  negative  correlations, 
and  further  multivariate  extensions. 

Finally  we  note  that  it  is  possible  to  introduce  non- 
stationarity  and  dependence  on  concomitant  variables  into  the 
sequence  by  multiplying  Xn  by,  say, 

r 

\ (n)  * exp{  l a . z . (n) } , 
j=0  3 D 

as  was  done  in  Cox  and  Lewis  (1966,  Ch.  3,  ii) . It  would 

be  of  interest  to  see  how  the  methods  given  in  Cox  and  Lewis 

extend  to  the  case  where  there  is  correlation  present  between 

the  X ' s . 
n 


31 


REFERENCES 


Bernier,  J.  (1970).  Inventaire  des  modeles  et  processus  stochastique 
applicables  de  la  description  des  deluts  journaliers  des  riviers. 
Revue  de  l'Institut  International  de  Statistique,  38,  50-71. 

Chemick,  M.  R.  (1977).  A limit  theorem  for  the  maximum  of  an 
exponential  autoregressive  process.  Tech.  Report  No.  14,  SIMS, 

Dept,  of  Statistics,  Stanford  University. 

Cox,  D.R.  (1962).  Renewal  Theory.  Methuen,  London. 

Cox,  D.R.  (1955).  Some  statistical  methods  connected  with  series 
of  events.  J.  Roy.  Statist.  Soc.  B,  1J,  129-164. 

Cox,  D.R.  and  Lewis,  P.A.W.  (1966).  The  Statistical  Analysis  of 
Series  of  Events,  Methuen,  London;  Wiley,  New  York. 

Downton,  F.  (1970).  Bivariate  exponential  distributions  of 
reliability  theory.  J.  R.  Statist.  Soc.  B,  32,  408-417. 

Feller,  W.  (1971).  "An  Introduction  to  Probability  Theory  and 
Its  Applications. " Wiley;  New  York. 

Gaver,  D.P.  (1972).  Point  process  problems  in  reliability. 

In  Stochastic  Point  Processes,  P.A.W.  Lewis  (ed. ) , Wiley,  New  York, 
775-806, 

Hammersley,  J.  and  Handscomb,  D.C.  (1964).  Monte  Carlo  Methods, 
Wiley,  New  York. 

Hawkes,  A.G.  (1972).  Mutually  exciting  point  processes.  In 
Stochastic  Point  Processes,  P.A.W.  Lewis  (ed.) , Wiley,  New  York, 
261-271'. 

Jacobs,  P . A.  (1978) . A closed  cyclic  queueing  network  with 
dependent  exponential  service  times.  J.  Appl.  Prob.  (to  appear). 

Jacobs,  P.A.  and  Lewis,  P.A.W.  (1977).  A mixed  autoregressive- 
moving  average  exponential  sequence  and  point  process  (EARMA  1,1). 
Adv.  Appl.  Prob.  9_,  87-104. 

Lawrance  A.J.  and  Lewis,  P.A.W.  (1977)  . A moving  average  exponential 
point  process  (EMAl) . J.  Appl.  Prob.  1£,  98-113. 

Lawrance,  A.J.  and  Lewis,  P.A.W.  (1978).  An  exponential  auto- 
regressive-moving average  process  EARMA (p,q) : Definition  and 
correlational  properties.  Naval  Postgraduate  School  Report 
NPS55-78-1 . 


32 


Lewis,  P.A.W.  and  Shedler,  G.S.  (1978).  Analysis  and  modelling 
of  point  processes  in  computer  systems.  Bull.  ISI  (to  appear). 

Neyman,  J.  and  Scott,  E.  (1972).  Processes  of  clustering  and 
applications.  In  Stochastic  Point  Processes,  P.A.W.  Lewis  (ed.), 
Wiley,  New  York,  646-681. 

Thorin,  0.  (1977) , On  the  infinite  divisibility  of  the  Pareto 

distribution.  Scand.  Actuarial  J.,  31-40. 

Wold,  H.  (1948) . Sur  les  processus  stationnaires  ponctuels. 
Collogues  Int.  Cent.  Natn.  Rech.  Scient.  1_3,  75-86. 


i 


33 


DISTRIBUTION  LIST 


No.  of  Copies 
1 


STATISTICS  AND  PROBABILITY  PROGRAM 
CPF  I CE  OF  NAVAL  RES  EARCF 
COLT  A 26 
ARLINGTON 

VA  22217 


CFFICE  CF  NAVAL  RESEARCH  1 

NHR  YORK  AREA  OFFICE 

715  BROADWAY  - 5 T F FLOOR 

ATTN:  CR.  ROEEP  GRAFT CN 

NEW  YORK,  NY  10033 


DIRECTOR 

OFFICE  OF  NAVAL  RESEARCH  BRANCH  OFF 
536  SOUTH  CLARK  STREET 
ATTN:  DEPUTY  A NO  CHIEF  SCIENTIST 

CHICAGO,  IL  60605 


1 


LI ERARY 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO 
CA 


92152 


1 


NAVY  LIBRARY 

NATIONAL  SPACE  TECHNOLOGY  LAB 
ATTN:  NAVY  LIBRARIAN 

BAY  ST.  LCL'IS 

MS  29522 


1 


NAVAL  ELECTRONIC  SYSTEMS  COMMAND  .1 

NAVE LEX  22C 

NATIONAL  CENTER  NO.  1 

ARLI  NC-TON 

VA  20360 


DIRECTOR  NAVAL  REAEARCF  LABORATORY 
ATTN:  LI  ERARY  ( ONRL ) 

CCCF  2 02  E 
WASHINGTON,  C.C. 

20375 


DEFENSE  DOCUMENTATION  CENTER  1 

CAFEKON  STATION 

ALEXANDRIA 
VIRGINIA  223 1 A 


TECtINI  CAl  INFORMATION  DIVISION  1 

NAVAL  RESEARCH  LABORATORY 

WASHINGTON,  D.  C. 


20375 


OFFICE  CF  NA  V A L RESEARCH 

SAN  FRANC  1SCO  AREA  OFFICE 

760  MARKET  STREET 

SAN  FRANC  ISC  C CALIFORNIA  94102 


TECHNICAL  LIERARY 
NAVAL  ORDNANCE  STATION 

INDIAN  HEAD  MARYLAND  2C64C 


NAVAL  SHIP  ENGINEERING  CENTER 
PHILADELPHIA 

DIVISION  TECHNICAL  LIBRARY 
PHILADELPHIA  PENNSYLVANIA  19112 


BIREAU  OF  NAVAL  PPSS0NN5L 
DEPART  PENT  CF  THE  NAVY 
TECHNICAL  LIERARY 
WASHINGTON  C.  C.  2 0370 


LIERARY  CCCE  0212 

NAVAL  PCSTGKACLATE  SCHCCL 

NCNTEREY  CALIFORNIA  9294C 


PRCF.  M . AECEL-HAPEED 
DEPART  HE  NT  CF  MATHEMATICS 
LNIVERS1TY  OF  NORTH  CAROLINA 
CHARLOTTE 
NC 

28223 


PROF.  T.  W « ANDERSON 
DEPARTMENT  OF  STATISTICS 
STANFORD  UNIVERSITY 

STANFORD  f CALIFCRNIA  943C5 


FRCF.  F.  J.  AN S COM BE 
DEPARTMENT  OF  STATISTICS 
YALE  UNIVERSITY 
NEK  HAVEN 

CCKNcCTICLT  C6520 


PROF.  L.  A.  ARCIAN 
INS  I T I III  S nr  IN1USTRIAL 
ALPINIST  RATION 
UNION  UJLl  l it 
SC.  I l Nit  I AD  \ , 

NEK  YORK  li.  300 


No.  of  Copies 

FRCF.  C.  R.  RAKFR  1 

DE  PAR  1 ME  NT  Or  STATISTICS 
UNIVERSITY  C.F  NCTRH  C A RCL  J NA 
CHAPEL  HILL, 

NCRTH  CAROLINA 
27514 


FRCF.  R.  E.  BECFFOFER 
CEFARTMEKT  CF  CPERATICNS  RESEARCH 
CCFNELL  UNIVERSITY 
IT FAC A 

NEW  YORK  14050 


FRCF  . N . J.  REPSHAC 
SCHOOL  OF  ENGINEERING 
UNIVERSITY  CF  CALIFORNIA 
IRVINE  , 

CA  LI  FCRNIA 
92664 


P.  J.  BICKEL 

DEPARTMENT  OF  STATISTICS 
UNIVERSITY  CF  CALIFCRN1  A 

BERKELEY  , CALIFCRNIA 

9472  0 


FFOF.  F.  W.  BLOCK 
D5  PART  ME  NT  CF  M A THE  MAT  I CS 
UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH 

FA 

1526 0 


PROF.  JOSEPH  BLUM 
DEPT.  OF  MATHEMATICS,  STATISTICS 
AND  COMPUTER  SCIENCE 
THE  AMERICAN  UNIVERSITY 
WASHINGTON 

CC  20016 


PROF  . R . A . BRADLEY 
DEPARTMENT  OF  STATISTICS 
FLORIDA  STATE  UNIVERSITY 

TALLAHASSEE  , FLORIDA  323C6 


PROF.  R.  E.  BARLOW  1 

OPERATIONS  RESEARCH  CENTER 
COLLEGE  CF  F NO  I N.:C  R IN  U 
UNIVERSITY  CF  CALIFORNIA 
BERKLEY 

CALIFCRNIA  94720 


MR.  C.  N FERN ETT  1 

NAVAL  COASTAL  SYSTEMS  LACCRATORY 

r.rur  p/ci 

FA  MAMA  CITY, 

FLORIDA 

22401 


No.  of  Copies 


PRCF . 1.  N.  THAT 

CGMPUTER  SCIENCE  / OPERATIONS 

research  CENTER 

SOUTHERN  METHL'L'IST  UNIVERSITY 

CALL  AS  , 

TEXAS  7 52  75 


FR  CF  * W.  F.  BUSCHKE 

DEFT.  CF  QUANTITATIVE 

BUSINESS  ANALYSIS 

UNIVERSITY  CF  SOUTHERN  CALIFORNIA 

LOS  ANGELES.  CALIFORNIA 

90007 


OR.  DERR  ILL  J.  BORDELON 

NAVAL  UNCEFWATER  SYSTEMS  CENTER 

CCJCc  21 

NEWPORT 

PI 

02840 


J.  E.  ECYER  JR 

DEPT.  OF  STATISTICS 

SOUTHERN  METHODIST  UNIVERSITY 

DALLAS 

TX 

75275 


DR.  J.  CHANCRA 
U.  S.  ARMY  RESEARCH 
F.  C.  C-OX  12211 
RESEARCH  TRIANGLE  PARK  , 
NGFTH  CAROLINA 
277C  t 


FRCF . H.  CHFRNCFF 

DEPT.  CF  MATHEMATICS 

MASS  INSTITUTE  CF  TECHNOLOGY 

CAMBRIDGE, 

MASSACHUSETTS  02139 


FFOF.  C.  GERMAN 

DE  FART  ME  NT  CF  CIVIL  ENGINEERING 
ANU  ENGINEERING  MECHANICS 
COLUMBIA  UNIVERSITY 
NEK  Y C f> K 

NEK  YORK  10027 


PRCF.  R.  L.  DISNEY 
VIRGINIA  POLYTECHNIC  INSTITUTE 
AND  STATE  UNIVERSITY 
DEFT.  CF  I K i:  U c T i\  J A I.  ENGINEERING 
ANU  OPERA  TICKS  RESEARCH 
BLACKS  BURG  , VA  24061 


MR.  J.  DCW11NC 

Dcli-N^E  LiH.  1 ST  ICS  STUDIES 

TNI  Ui'.FA TICK  LXCHAM'.f- 

ARMY  lUGlSTUS  MANAGEMENT  CENTER 

FUKT  l:»;  » 

VIRGINIA  20390 


1 


1 


1 


1 


1 


1 


1 


1 


1 


No.  of  Copies 
1 


FRCF  J.  C.  FSARY 

DEFT.  OF  OPERATIONS  RESEARCH  AND 
AUMIN'ISTRAT  IVE  SCIENCES 
NAVAL  PO  £ 1 GRACIA  7 i:  SCHGCL 
MONTEREY  , 

CALIFORNIA  9394D 


CF.  V . J.  FISCHER 
DEFENSE  COMMUNICATIONS  AGENCY 
I860  WIEPLE  AVENUE 
REST  CN 
VIRGINIA 
22C73 


FRCF.  D„  P.  GAVER 

CEPT.  CF  OPERATIONS  RESEARCH 

NAVAL  POSTGRACLATfc  SChCCL 

MONTEREY 

CA 

93940 


MR.  GENE  H.  GLEISSNER 
AFFLIEC  MATH  EM AT ICS  LABCRATCRY 
CAVID  TAYLCR  NAVAL  SHIP  RESEARCH 
ANT  DEVELOP PENT  CENTER 
EETHcSDA 

MD  20084 


FROF.  S.  S.  GCPTA. 
DEPARTMENT  Of  STATISTICS 
PLRCU0  UNIVERSITY 
LAFAYE  T T E 
INCIANA  47907 


PROF.  C.  1.  HANSON 
DEPT  OF  NATH.  SCIENCES 
STATE  UNIVERSITY  OF  NEW  YORK, 
BINGHAMTON 
BINGHAMTON 

NY  13901 


PROF.  F.  J.  HARRIS 

CEPT.  CF  ELECTRICAL  ENGINEERING 

SAN  D! EGG  STATE  UNIVERSITY 

SAN  DIEGC 

CA 

92182 


PROF.  L.  P.  HER BACH 
DEFT.  CF  OPERATIONS  RESEARCH  AND 
SYSTEMS  ANALYSIS 

POLYTECHNIC  INSTITUTE  OF  NEW  YORK 
ORCQKLYN 

NY  11201 


FRCF.  M.  J.  HMCH 
DEPARTMENT  CF  ECONOMICS 
VIRGINIA  PCLYTfCl-MC  INSTITUTE 
AND  STATE  IMVEHS1TY 
B LAC  RS'UJRG  . 

VUU.1MA  2A061 


1 


1 


1 


1 


1 


1 


1 


1 


No.  of  Copies 

1 


PPOF.  W . P.  H IRSCH 
INSTITUTE  Cl-  RATHE  V ATI  C A L SCIENCES 
NE  R YORK  UNIVERSITY 
Mi R YORK 
NEU  YORK  1C453 


FRCF.  0.  L.  I GLEHART 
DEPARTMENT  CF  OPERATIONS  RESEARCH 
STANFORD  UNIVERSITY 
STANFORD  , 

CAL  1 FORM  I A 
54  -SO 


FRCF  . J.  e.  K A CANE 
DEPARTMENT  CF  STATISTICS 
CAPNHG  IG-MELLGN 
FITTS  EURCH  , 

PE  NNSYLVAMA 
15213 


DR.  RICHARD  LAU 
Cl FCCTOR 

CFTICE  CF  NAVAL  RESEARCH  ERANCH  OFF 

1 C30  EAST  GREEN  STREET 

PASADENA 

CA  $1101 


DF.  A.  R.  LAUFER 
DIRECTOR 

OFFICE  OF  NAVAL  RESEARCH  BRANCH  OFF 

1030  EAST  GREEN  STREET 

FASACENA 

CA  $1101 


PROF.  N.  LEADBETTER 
DEPARTMENT  OF  STATISTICS 
UNIVERSITY  OF  NORTH  CAROLINA 
CHAPEL  HILL 
NOFTH  CAROLINA  27514 


CR.  J.  S.  LEE 

J.  S.  LEE  ASSOCIATES,  INC. 

20C1  JSFFEBSCN  DAVIS  HIGHWAY 

SUITE  80  2 

ARLINGTON 

VA  22202 


FRCF.  1.  C.  LEE 
DEPARTMENT  CF  STATISTICS 
VIRGINIA  POLYTECHNIC  INSTITUTE 
AND  STATE  UNIVERSITY 
BLACKSBURG 

VA  24061 

FFCF.  R.  S.  L FV  GNK'ClRTli 
CEI'T.  CF  INDUSTRIAL  ANC  SYSTEMS 
LNCH'P  F.K  TNG 
LNIV'NSnV  CF  FLCFIOA 
UAINSV1LI.E  , 

FL  Cl>  1 DA  32(11 


1 


1 


1 


1 


1 


1 


j 


I 


FRCF.  P.  A.  to . LEWIS 

CEPT.  CF  OPERATIONS  RESEARCH  AND 

ADMINISTRATIVE  SCIENCES 

NAVAL  FCST  CRUCIATE  SCFCOL 

MGNTEREY  , CALIFORNIA 

93SA0 


FRCF  G.  LI EEcRMAN 

STANFORD  1N1VHRSITY 

DEPARTMENT  CF  OPERATIONS  RESEARCH 

STANFORD  CALIFORNIA  S4305 


DR.  JAPES  R.  MAAR 
NATIONAL  SECURITY  AGENCY 
FORT  MEADE  , MARYLAND 
20755 


FRCF.  R.  to.  MACSEN 
DEPARTMENT  CF  STATISTICS 
UNIVERSITY  OF  MISSOURI 
COLUMBIA 
MO 

65201 


DR.  N.  R.  MANN 
SCIENCE  CENTER 

ROCKWELL  INTERNATIONAL  CGRFCRAT I CN 
F.C.  BOX  108  £ 

TFCUSAND  C/-KS  » 

CALIFORNIA  913CC 


CR.  K.  h.  MARLOW 

PROGRAM  IN  LOGISTICS 

TFE  GEORGE  toA  SHI  NOT  ON  UNIVERSITY 

70?  22NC  STREET  , N.  W. 

WASHINGTON  , 0.  C. 

2003  7 


PROF.  E.  MASRY 
OEFf.  APPLIED  PHYSICS  AND 
INFORMATION  SERVICE 
UNIVERSITY  OF  CALIFORNIA 
LA  JCLIA 

CALIFORNIA  92093 


DR.  PRICE  J.  MCCONALD 
SCIENTIFIC  DIRECTOR 
SCIENTIFIC  LIAISON  CROUP 
OFFICE  CF  NAVAL  RE  Sr Al<  CF 
AMERICAN  EMBASSY  - TOKYO 
AFC  SAN  f R /NCI  £ CO  96503 


PROF.  J . A.  MLCKSTADT 
DEFT.  CF  OFFRAIKNS  RESEARCH 
CCFNI:!I.  INIVERSITY 
IT  FAC.  A , 

NOW  YLFK 


No.  of  Copies 


CR  . JANET  M . M VMRE 
THE  INST17L7E  OF  DECISION  SCIENCI 
FOR  BUSINESS  ANC  PUBLIC  POLICY 

Claremont  men's  college 

CLAREMONT 

C A 51711 


MR  . P . NISSELSCN 
BIPEAL  GF  TEL:  CENSUS 
ROOM  2C2  5 

FPCERAL  f'JILCII'G  3 
NASH1  NG7CN  , 

0.  C.  2 03- 


MISS  E.  S.  ORLEANS 
NAVAL  SE/  SYSTEMS  COMMAND 
(SEA  03  F ) 

RM  IOSCO 

ARLINGTON  VIRGINIA  20360 


FRCF.  C.  E OWEN 

DEPARTNENT  CF  STATISTICS 

SOUTHERN  METHOCIST  UNIVERSITY 

CALLAS 

TEXAS 

75222 


PROF.  E.  PAR ZEN 
STATISTICAL  SCIENCE  DIVISION 
STATE  LNIVERSI7Y  OF  NEW  YORK 
AT  BUFFALO 
AMFERST 

NEW  YORK  14226 


DR.  A.  FETRASOVITS 

ROOM  20  7 E , FOOD  AND  CRcG  ELDG. 

TUN’Nc  Y • S PASTLRc 

CTT01-.A  * ONTARIO  K1A-CL2  » 

CANADA 


FRCF.  S.  L.  PFCENIX 
SIELEY  SCHOOL  OF  MECHANICAL  AND 
AEROSPACE  ENGINEERING 
CORNELL  UNIVERSITY 
ITHACA 

NY  14850 


DP.  A.  L W POWELL 
Cl  RECTOR 

OFFICE  OF  NAVAL  RESEARCH  BRANCH  OFF 

4S5  SUMMER  STREET 

BCSTCN 

MA  02210 


MR.  F.  R.  FPJCFI 

CODE  P2'i  OPE  RAT  ION  SL  TEST  AND  DNRS 
EVALUATION  FilKCE  (UPTCVFDiU 
NC NT  IN  K , 

V IRC.  1 NIA 

20;.ou 


PROF.  M.  L.  PURI 
DEFT.  CF  MATH2M AT  I CS 
P.G.  BOX  F 

INDIANA  LN  I VERS  IT  Y FOUNDATION 
E l CO  MI NGTCN 

IN  47401 


FROF.  H RCEBINS 
DEPARTMENT  Of  MATF2MAT  ICS 
CCLUI'EIA  UMVEPSITY 
NEW  YORK, 

NEW  YORK  1CJ27 


PFOF.  M R0SSN3 1.ATT 
DEPARTMENT  CF  MATHEMATICS 
UNIVERSITY  CF  CALIFORNIA  SAN  DIF.GO 
LA  JCLLA 
CALI FCRNI  A 

92093 


PROF.  S.  M.  PCSS 
COLLEGE  CF  ENGINEERING 
UNIVERSITY  CF  CALIFORNIA 
BERKELEY 
CA 

94720 


PROF.  I RUBIN 

SCHOOL  OF  ENGINEERING  AND  APPLIED 
SC  IENCE 

UNIVERSITY  OF  CALIFORNIA 
LCS  ANGELES  , 

CALIFORNIA  90024 


FPCF.  I.  R.  SAVAGE 
CEPARTMENTOF  STATISTICS 
YALE  UNIVERSITY 
NEW  HAVEN, 

CONNECTICUT 

C6LS2U 


FPCF.  L.  L.  , SCHARF  JF 
DEPARTMENT  CF  ELECRICAL  ENGINEERING 
COLOR  A CO  STATE  UNIVERSITY 
FT  . CCLLINS  , 

CCLORACO 

60521 


PROF  . R.  SFRFLING 
DEPARTMENT  CF  STATISTICS 
FLCRIOA  STATE  UNIVERSITY 

TALLAHASSEE  FLORIDA  32306 


PROF.  I,  . R.  SCMLCANY 
DEPARTMENT  OF  STATISTICS 
SU1.TH&RK  ME  THCCI  ST  UNIVERSITY 


No.  of  Copies 


PROF.  C.  C.  SISGMUND  * 1 

CEPT  . OF  STATIST  ICS 
STANPCRC  UNIVERSITY 
STANFORD 

CA  54305 


FRCF.  M.  L.  S FOQMAN  , 1 

DEPT.  CF  ELECTRICAL  ENGINEERING 
POLYTECHNIC  INSTITUTE  CF  NEW  YORK 
BRCCKLYN* 

NEW  YORK 
11201 


PROF.  N.  SING  POP.  WALLA 
CEPT.  CF  OPERATIONS  RESEARCH 
THE  GEORGE  WASHINGTON  UNIVERSITY 
7C7  2 2 NO  ST  . N.  W. 

WASHINGTON.  C.  C.  _ , 

20052 


DR.  A.  L.  SLAFKCSKY  1 

SCIENTIFIC  ADVISOR 

COMMANDANT  OF  T F c MARINE  CORPS 

WASHINGTON  , 

D.  C. 

20380 


PR.  CHARLES  S.  SMITH  1 

C AS  C ( I £ L ) * 

PENTAGON 

WASHINGTON 

CC 

20301 


CP.  C.  E.  SMITH  1 

DESMA7ICS  INC. 

P.C.  BOX  618 
STATE  COLLEGE 
PENNSYLVANIA 
16C01 


PROF.  W.  L.  SMITH 
DEPARTMENT  OF  STATISTICS 
LNIVFRSITY  CF  NORTH  CAROLINA 
CHAPEL  HILL 
NORTH  CAROLINA  27514 


FRCF.  H SCLCKCN 
CEPARTMENT  CF  STATISTICS 
STANHCRG  UNIVERS  ITY 
STANFORD  » 

CALIFORNIA 

54305 


1 


MR.  GLENN  F.  STAHLY  1 

NATIONAL  SECURITY  AGENCY 
FORI  Mi  FACE  . 

MARYLAND  20755 


No.  of  Copies 
1 


NR.  CAVIC  A.  SWICK 
ADVANCED  PROJECTS  GROUP 
CODE  8 1 C2 

NAVAL  RESEARCH  LAB. 

NASHINGTCN 

CC 


20375 


NR.  WENDELL  G.  SYKES 
ARTHUR  C.  LITTLE,  INC. 
AC CRN  PARK 
CAR-BRIDGE 
HA 


02140 


1 


PROF.  J.  P.  THCMPSON  „ „ 1 

DEPARTMENT  OF  MATHEMATICAL  SCIENCE 
RICE  UNIVERSITY 
H COS  TON  , 

TEXAS 

77001  1 


PROF.  W.  A.  TI-CNFSCN  1 

DEPARTMENT  OF  STATISTICS 

UNIVERSITY  CF  NISSGUR1  1 

COLUMBIA  , 

MISSOURI 

65201 


FRCF.  F . A . T ILLMAN  1 

DEPT.  CF  INDUSTRIAL  ENGINEERING 

KANSAS  STATE  UNIVERSITY 

MANHATTAN 

KS 

66506 


PRCF  J.  V.  TU KEY  1 

DEPARTMENT  CF  STATISTICS 
PRINCETON  UNIVERSITY 

PRINCETON  , N.  J.  08540 


PROF.  A . F . VSINOTT 

DEPARTMENT  CF  OPERATIONS  RESEARCH 

STANFORD  UMVESST1TY 

ST ANFCPC 

CALIFORNIA 

94305 


DANIEL  H.  WAGNER 

STATION  SUL  ARE  ONE  1 

FACLI  , PENNSYLVANIA 

15301 


PRCF.  GRACE  WAHPA 
run . CF  STATISTICS 
LN  I Vf.P.  S 1TY  CF  WISCONSIN 
MADISON 

W1 


1 


53706 


No.  of  Copies 
1 


PRfF.  K.  T.  WALLEMUS 

DE  FAR  THEN  7 OF  MAThEMAT  ICA  L SCIENCES. 

CL  EM SON  UNIVERSITY 

Clfc'MSCNt 

SOUTH  CARCLINA  29631 


PPCF . G.  S.  WATSCN 
DE  FART ME NT  OF  STATISTICS 
PRINCETON  . N.  J.  C854G 


1 


PROF.  BERNARD  hlDRCW 
STANFCRC  ELECTRGNICS  LAE 
STANFORD  UNIVERSITY 
STANFORD 
CA 


94305 


FRCF.  G.  f.  W F ITEFCUSE 

DEPT.  OF  IN  CL  SIR  IAL  ENGINEERING 

LEFIGF  UNIVERSITY 

BETHLEHEM 

PA 

18015 


FFOF.  S.  ZACKS 

DEPT.  O'-  MATHEMATICS  AND  STATISTICS 
CASE  WESTERN  RESERVE  UNIVERSITY 
CLEVELAND. » 

CHIO 

44IC6 


PFOF.  M.  Z IA-HAS SAN 
DEPARTMENT  CF  INDUSTRIAL  AND 
SYSTEMS  ENGINEERING 
ILLINOIS  INSTITUTE  CF  TECHNOLOGY 
CUCAGC 

IL  60616 


HEAD,  MATH.  SCI  SECTION 
NAT.  SCIENCE  FOUNDATION 
18C0  G STREET,  N.W. 

WASHINGTON,  C.C. 

205  50 


PROF.  A.  F.  ANCRUS 
CEPT.  CF  CP 

NAVAL  FCSTORACLATE  SCHCCL 

MONTEREY 

CALIFORNIA 

53940 


FFOF.  C.  R . 13ARR 
Pfc  FT  . CF  CP 

NAVAl  POSTGRACLATE  SCHOCL 
MON TURLY 
CAL  It  CNN  I A 


1 


1 


1 


1 


1 


1 


1 


93940 


No.  of  Copies 
1 


PPCF . PETER  BICCMFIELC 
STATISTICAL  CcPT. 

PRINCETON  UNIVERSITY 
PRINCETON,  N.  J. 

C8540 


PPCF.  G.  G.  BROWN 
CE  FT  . CF  OP 

NAVAL  FQSTGP ACUATE  SCFCCL 

MONTEREY 

CALI FCPNIA 

S3  940 


R.  W.  PUTTERWCRTH 
SYSTEMS  EXPLORATION 
WEBSTER  ST. 

MONTEREY 

CALIFORNIA 

939  AO 


DR.  JAMES  CAPRA 
7218  CELFIELC  STREET 
CHEVY  CHASE 
MAPYLANC 

20015 


CR.  C.R.  COX 
DEPT.  CF  MATHEMATICS 
IMPERIAL  COLLEGE 
LONDON  SW7 

E NG  LA  N 


1 

DEFENSE  OGCUMENTATI CN  CTR. 

CAMEPCN  STATION 

ALEXANDRIA 

VIRGINIA 

22314 


ECCN.  AND  NAN.  SCI.  CTR. 
NORTHWESTERN  UNIV. 

EVANSTON 

ILLINOIS 

602  Cl 


MAN.  SCE.  RES.  CTR. 

FACULTY  OF  CUM.  AND  BUS.  ADMIN. 
LNIV.  CF  EPIIISH  COLUMBIA 
VANCOUVER 

BRITISH  CCLUMEIA  V6T  1W5  CANADA 


CR.  M.  CN  ASS 
MATN.  DL  FT  . 

NORTHWESTERN  UNIV. 

EVANSTON 

ILLINOIS 

602  01 


No.  of  Copies 

1 


CR.  R . BLASHCFP 
BIOMAT  HE  VAT  ICS 
UMV.  CP  C /LI  F . 
LOS  ANGELES 
CALIFORNIA 


90024 


PROP.  OECRGE  S.  FISHMAN 
LMV.  CF  NCRTI-  CAROLINA 
CUR.  IN  OR  AND  SYS.  ANALYSIS 
PHILLIPS  ANNEX 
CHAPEL  HILL*  NORTH  CAROLINA 

20742 


DR.  R.  GNA PACES  I KAN 
CELL  TELEPHONE  LAB 
HO  LI-'DE  L t N • J. 


07733 


DR.  A.  J.  GCL  OMAN 
CH I E F C R 

DIV.  2C5.  C2 1 ADMIN.  AA26 
US.  DEPT.  CF  COMMERCE 
WASHI ngton ♦ c.c. 

20234 


CR.  H.  PIGGINS 

53  BONN  I*  POST FACH  5flS 

NASSeSTRASSE  2 


WEST  GERMANY 


DP.  P.  7.  HOLMES 
DEPT  . OF  MATH. 
CLEMSCN  LNIV. 
CLEM  SON 
SCUTH  CAROLINA 


29631 


CR  . J . A . FOCKI 
EELL  TELEPHONE  LAeS 
HO  LMDS L 
NEW  JERSEY 


077  33 


CR.  ROBERT  IWCKS 
MATH.  DEPT. 

WEST  JNCHQUSE  RES  . LABS 
CHURCHILL  130  FO 
PITTSBURGH,  PENNSYLVANIA 

152.15 


CR.  0.  L.  I CL  EH ART 
CEPT.  CF  C.H. 
stand  run  LMV. 

STAN) UNO 
CALI  l;L I N 1 A 


1 


1 


1 


1 


1 


1 


1 


1 


943  05 


No.  of  Copies 


CP.  PATRICIA  JAC CBS 

CR  DEPT. 

NAVAL  POST GRADLA7E  SCHCCL 
MONTEREY 

CALIFORNIA 

93940 

1 

DR.  H.  KOEAYASFI 

IB  M 

YCRKTCl’N  HEIGHTS 

NEW  YORK 

10598 

1 

CR.  JOHN  L EFOCZKY 

STATISTICS  CtPARTFENT 

CARNEG  Iifl-H’E  LLCN  UNIVERSITY 
FITTS  BURGH 

PENNSYLVANIA 

15213 

1 

L IERARY 

CC C ^ 5 5 

NAVAL  PCSTGRACUATE  SCHCCL 
MONTEREY 

CALIFORNIA 

93940 

1 

DR.  A.  L ENGINE 

1020  GUINEA  ST. 

FALQ  ALTC, 

CALI  FORN  1A 

943  01 

1 

CR.  J.  MACCUSEN 

UNIV.  CF  CALIF. 

LOS  ANGELES 

CALIFORNIA 

90024 

1 

FPCF.  K.  T.  MARSHALL 

CELT.  CF  CF 

NAVAL  POSTCRACUATE  SCHCCL 
MONTEREY 

CALI FCPNIA 

93  9 4 0 

* 

1 

DR.  M.  MA2UMCAR 

MATH.  DEPT. 

ESTiNCFOLSi  RES.  LABS 
CHURCHILL  3CFC 

FITT  SUIRC-H 

PENNSYLVANIA 

15235 

1 

DR.  t. F ON  F.  MCGINNIS  1 

SCHOLL  Cf  INC.  A NC  SYS.  ENG. 

GEORGIA  INST.  OF  TECH. 

ATLANTA 
CLCRuI A 


30332 


T 


t 


CR.  D.  R.  MCNEIL 
DEFT.  CF  STATISTICS 
PR  INC E TON  UNIV. 
PRINCETON 
NE In  JERSEY 


08540 


PPCF . P.  R.  MILCF 
CC FT  , CF  OR 

NAVAL  FOSTORACUATE  SCFCCL 

MONTEREY 

CALIFORNIA 


93940 


CR.  F.  MOSTELLER 
STAT.  CcFT. 

FARVARC  CMV. 

CAMBRICGE 
MASSAC FU SETTS 

02139 


FRCF . R.  R.  READ 
DEFT  . CF  CP 

NAVAL  FCSTGRACUATE  SCFOOL 

MONTEREY 

CALIFORNIA 


93940 


DR  . M . REI SER 
I EM 

THOMAS  J.  WATS  C N RES.  CTR. 
YCRKTOWN  HEIGHTS 
NEW  YCRK 


10598 


OEAN  CF  RESEARCH 
CODE  013 

NAVAL  FCSTGFACLATE  SCFCOL 
MONTEREY 

CALIFORNIA  93940 


FRCF.  F.  R.  RICHARDS 
DEFT.  CF  OR 

NAVAL  FCSTGFACLATE  SCFCOL 

MCNTcRY 

CALIFORNIA 

93940 


OR.  J.  RICRCAN 
CEFT  . OF  MATHEMATICS 
ROCKEFELLER  UNIV. 

NEW  YORK 
NEW  YORK 


DR.  LINUS  SCHR/GE 
LNIV.  C F CFICAGO 
CFAP.  SC  FOOL  OF  hi S • 
50  So  GREENWICI.  AVt. 
CHICAGO.  ILLINOIS 


.100  21 


No.  of  Copies 
1 


1 


1 


1 


1 


1 


1 


1 


1 


60637 


No.  of  Copies 


DR.  PAUL  SCHhEITZER 

TMCMAS  J.  WATSCN  RESEARCH  CTR9 


P C ST  OFFICE  BOX  218 
YGRKTONN  FLIGHTS 
NEW  YCFK 


CR.  RICFAPC  SCPENSCN 
CODE  303  NPRDC 
2?1  CATALINA  BLVD. 
SAN  C1EGC 
CALIFORNIA 


FFCF.  M.  G.  SOVEREIGN 
DE  FT  . CF  OP 

NAVAL  F0S7GRACUATE  SCFCCL 

MONTEREY 

CALI FCRNIA 


10598 

1 


92152 

1 


93  940 


CR.  V.  S FlNIV  AS  AN 

CR  AD  LA  Tc  SCHCCL  CF  BUSINESS 

STANFORD  UNIVERSITY 

STANFORD 

CALI FCRNIA 

94305 


DR  . R.  M.  STARK  1 

STATISTICS  ANC  CCNPUTEF  SCI. 

UNIV.  CF  DELAWARE 

NEWARK 

CELAWARE 

19711 

PFOF.  RICFARC  VANS  LYKc  1 

RES.  ANALYSIS  CCRP. 

BEEC  l-WCOC 

CLD  TAFFFN  FG/C 

GLEN  CCVLt  NEW  YORK 

11542 


PRCF • JOHN  W.  7UKEY  1 

FINE  FALL 
F P 1 KC-C TC1N  UNIV. 

PRINCETON 
NEW  JLRSEY 

08540 


CR  . TFCMAS  C.  VARLEY  1 

CFFICE  OF  NAVAL  RESEARCH 

CODE  4 24 

ARLINGTCN 

VA 

22217 


FRCF.  G . WATSCN 
F 1NE  II  h I 
FRlNCf  I.:N  UMV. 
PR  I NCI  II  N 
NEW  JERSEY 


1 


COS  40 


Dr.  Roy  Welsch 
M. I ,T. , Sloan  School 
Cambridge,  MA  02139 


No.  of  Copies 


1 


Dean  of  Research  012 

Naval  Postgraduate  School  1 

Monterey,  Ca.  93940 

Professor  D.  P.  Gaver  50 

Code  SSGv 

Naval  Postgraduate  School 
Monterey,  Ca.  93940 

Professor  Peter  Lewis 
Code  55 Lw 

Naval  Postgraduate  School 
Monterey,  Ca.  93940 

R.  J.  Stamp fel 
Code  55 

Naval  Postgraduate  School 
Monterey,  Ca.  93940 


